NASA TND-1540 



TECHNICAL NOTE 

D-1540 


REMARKS ON HILL’S LUNAR THEORY. 
PART I 


Karl Stumpff 

Goddard Space Flight Center 
Greenbelt, Maryland 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON March 1963 








REMARKS ON HILL’S LUNAR THEORY. 
PART I 
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Goddard Space Flight Center 
SUMMARY 

G. W. Hill, by neglecting solar parallax and the eccentricity of the 
solar orbit, presented a particularly simple form of the differential 
equation of the restricted three-body problem for the motion of a mass- 
less satellite. In the present paper, the Hill equation is modified to give 
a third order differential equation for r— the planetocentric distance of 
the satellite. And this equation can be solved by iteration (if the satel- 
lite orbit is considered to be a disturbed Kepler ellipse). This new 
solution is not only suitable for computer applications but it is inde- 
pendent of the coordinate system and is valid for both a fixed and a 
rotating system whereas Hill's solution was limited to simple-periodic 
orbits in a rotating system. The discussion of Hill's lunar problem 
provides a simple example of a method applicable also to more difficult 
problems such as the restricted three-body problem, in which conic- 
section orbits can be considered good approximate solutions. 
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REMARKS ON HILL’S LUNAR THEORY. 
PART I 


by 

Karl Stumpff* 

Goddard Space Flight Center 


INTRODUCTION 

By discounting the solar parallax, G. W. Hill 1 ^ gave a particularly simple form to the differential 
equation of the restricted three-body problem for the motion of a massless satellite around a planet 
in a circular orbit about the sun. This relatively simple problem offers tempting considerations and 
experiments, which will be discussed here. In this paper, it will be shown that the Hill Equation can 
be modified so that a differential equation of the third order for r —the planetocentric distance of the 
satellite— results, in which (besides r and its derivations through the third order) only the Jacobian 
constant occurs. And this equation can be solved by iteration if the satellite orbit is considered to 
be a disturbed Kepler ellipse. 


HILL'S DIFFERENTIAL EQUATIONS 

Hill's differential equations can be written as: 


where 


= -*(£-$) ■ 
’i + = "77 ’ 


r 2 - £2 + Tj 2 , 




V 


> (1) 



Figure 1 —Geometric representation 

Here <f, v are the right-angle planetocentric coordinates of the 0 f a mass | e ss satellite sofa planet 

assumed massless satellite S in a coordinate system (Figure 1) P- 


•NAS-NASA Research^Associate; Professor Emeritus, Gottingen University. 

tHill, G. W "Researches in the Lunar Theory,” in The Collected Works of George William Hill, Vol. t, Washington, D. C.: Carnegie 
Institution, 1905 • 


1 



which revolves with the constant angular velocity n = 1. And we have 


n 



= 1 


( 2 ) 


for the mean motion of the planet P in its circular orbit around the sun S, which is to be assumed at 
a great distance, a, along the negative axis* 


The units of mass are to be so selected that Equation 2 is valid. This is the case when we set 
k = 1 and a 3 = nij^ + m where m x is the mass of the sun and m the mass of the planet. In addition, the 
unit of mass is selected so that m = 1. In the case of the motion of a satellite, or moon, about the 
earth, £ 330,000 and the unit of time is equal to 1A = 58. d 13244. It also follows that 
a = + m £ 69. Thus Equation 1 is obtained if, in the differential equation of the restricted three- 

body problem, m t and a are allowed to go towards infinity so far that mj/a 3 - 1, but 
rr^/a 5 , m t /a 7 , • •• - 0. This simplification of the problem means that the sun, as a disturbing body, 
is removed to any great distance desired (so that its parallax at the earth-moon distance becomes 
unnoticeably small), but that its disturbing effect upon the satellite’s motion (which is proportional 
mja" 3 ) about the earth remains unchanged. 


Hill's equations (Equation 1) cannot simply be transformed into the equations of the undisturbed 
motion of the satellite, since the mass of the disturbing body (the sun) does not occur explicitly in 
them, but rather is implicitly contained in the definition of the unit of length. The equations can, 
however, be given a form which makes it possible to combine the disturbed and undisturbed motion. 
If we write 


where 


* - 2 V = *(l 2k) , • 

V + 2£ = 7,(l -± - k) , 


(3) 


> 


£ 2 + 7] 2 - r 2 


J 


then these equations assume the form of Equation 1 for k = 1, whereas letting k = 0 yields the equa- 
tions of undisturbed motion (in a coordinate system revolving with the angular velocity n = 1). 

Equations 3 have the Jacobian integral 

i 2 + V 2 = r 2 ( 1 -k) + y+ 3k,f 2 - C . (4) 
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THE DIFFERENTIAL EQUATION FOR r(t) 

The two differential equations of the second order (Equation 3) form a fourth-order system. In 
place of them, a single fourth-order differential equation can be written for a suitable variable — per- 
haps r = V<f 2 + tj j . With the help of the Jacobian integral, the order of this equation can be decreased 
by one, so that we have a relation of the form 

f (lr # r , f , r ; C ) = 0 . (5) 

This same consideration applies mutatis mutandis for the complicated probleme restreint , which we 
will discuss, but first we will make this experiment with the simple Hill problem. 

If we differentiate 


r 2 = £ 2 + T } 2 


( 6 ) 


three times, and if we eliminate 77 by using Equation 3 and £ 2 + ^by using Equation 4, we obtain 

rf = f + 7]7] , 


1 


r r + r 2 = 2(£ r\ - r \£ ) + 2r 2 ( 1 - k ) + -p + 6k - C 


rr’ + 3rr = 6k( 2£ f - <f 77 ) " “ 4k r r . 


> 


(7) 


Equations 6, 7 and 4 form a system 

<£ 0 = r 2 = £ 2 + v 2 » (8a) 

4> x - rv = f j ( 8 b) 

<*> 2 = r r + P "T “ 2r 2 (1 -k) + C = 2(£ 7 , - tj£) + 6 k £ 2 . ( 8 c) 

<£ = r r’ + 3r r + — r + 4kr f = 6 k( 2 £f-£ 7 ?) , ( 8 d) 

3 r * 

<p 4 = r 2 (l - k) + ~~C = <f s + V 3 - 3k £ 2 , ( 8 e) 


in which five functions <P 0 , ■■■ <t> A of r, r, r, F, and C appear on the left and appear on the right as 
functions of the four quantities #, v, <f, V- Therefore, to obtain the desired equation in the form of 
Equation 5, it suffices to eliminate the latter four quantities from Equation 8 . 
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ELIMINATION 


It follows from Equations 8b and 8c, and by taking Equation 8a, that 

= 2<j> l g - <p 2 7] + 6k£ 2 7) 

W 0 V = 2*^7? + <fi 2 £ - 6kf 3 . 

If we multiply Equation 9a by £, Equation 9b by and add, we get 

w 0 (€ 2 + 7 2 ) = + ff n - if) (<t> 2 - 6k f 2 ) 

or, since from Equation 8c it follows that <£ i) - -6kf 2 ) , 

^>o + V) = + ?(<^ 2 ~ 6k£ 2 ) 2 . 

If Equations 9 and 10 are substituted into Equations 8d and 8e, there results 

*0^3 = 6k [ 2 *i^ 2 ~ €v(<P 0 + <f> 2 ~ 6k£ 2 )] , 

+ l(b “ 6k^ 2 ) 2 - 3k<£ 0 £ 2 , 

and from Equation 8a 


(9a) 

(9b) 


(10) 


(11a) 

(lib) 


4> o e = + t^) 2 • 


(11c) 


The 6kf 2 and the 6k^ terms can easily be eliminated from Equations 11 by letting x = 6k£ 2 , y = 
for if we set 


<£ 0 + “ X = \ , 


from Equation lib we get 

4 ^<A - ^i 2 ) = (* - ^ 0 ) 2 - »o(* 0 ' + - M = k 2 + 4> 0 2 - 20 o (* 0 + * 2 ) . 

The above equation means that 


k 2 - 4>i? + 2<p Q [4> 2 + 24> a ) - 4^, 2 
is dependent only upon 4> 0 , 4> v 4> 2 , <t> 4 ; i.e., upon r, r, ?, c. 


(12) 


(13) 
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From Equations 11a, 11c, and 12 we have 


y\ - 2 4> 1 ( <t> 0 + <t> 2 ~ K ) ~ ^0^3 • 


y 2 - "(<£ 0 + <t> 2 ~ X) 2 + 6k 4> 0 (<t> 0 + <p 2 - , 


from which, we obtain a fourth-order equation in X by eliminating y: 


X 4 + cijX 3 + a 2 X 2 + a 3 X + a 4 - 0 


(14) 


where the coefficients of X are dependent only upon 4> 0 , • ■ • tf> + : 


a, = 2[<*> 0 (3k - 1) - <t> 2 ] , 


a 2 = 4 <p* - 9k 2 4 > 0 2 + -4 cl* , 


a 3 ~ 4 ^i [* 0*3 " 2< ^i K + ^2)] ’ 

% = [*•*. - + ^)] 2 • 


(15) 




It is worthy of note that 'f occurs only in <£ 3 , specifically in the combination ^ 0^3 - 2^1 K ^ ^ 2 ) * 
which, in turn, appears only in a 3 and a 4 . Therefore, if we set 


f - <V3 " ^lK + ^ 2 ) > 


(16) 


we can also write Equation 14 as 


f 2 + 4 <P^k( + X. 2 (a 2 + + k 2 ) - 0 . 


The above equation has the following solution: 


f = - 2<£, X ± X y^j 2 - (a 2 +a t X + X 2 ) 


or, if a. 1? a 2 are e?q?ressed by Equation 15, 


f - - 24> x k± k y^9k 2 0 o 2 “ [ X + <p 0 { 3k - 1 ) - <£ 2 ] 2 ■ 


( 17 ) 


Equation 17, together with Equations 13 and 16, provides the formal solution of the problem. 
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APPLICATION 


If we consider the motion of the satellite as a disturbed Kepler motion, it is practical to recall 
that for k = 0 the third of Equations 7 transforms into the relation 


r r’ + 3f r 



0 , 


(18) 


which can be integrated in a closed form. For, according to Equation 8d, the disturbed motion is 


r r + 3r r 



<p 3 - 4k 4> x . 


But from Equation 16, 


(* 3 - 4k ^i) = f'+ 20,(0 O + 0 2 - 2k 0 O ) ; 


therefore 


4 > q {t r + 3r r + -jfj) ~ 2 0, (0 O + 0 2 " X - 2k 0 O ) ± X ]/9k 2 0 O 2 -[<£„+ - X - 3k0 o ] 2 

The right side of this equation becomes zero for k = 0, since, according to Equation 12, 
0 O + 0 2 - \ = x = 6k£ 2 also contains k as a factor. 

For the case in which the satellite's orbit can, to first approximation, be considered a Kepler 
ellipse, the right side of Equation 19 is small, i.e., of the order of the perturbation. We then have, 
for k = 1: 


00 (r r + 3f r + ~J = 20, (0 2 - 0 O - x) ± X VW» 0 2 ~ {<t> 2 - 20 o - X) 2 


with 


l (20) 


k* = 


= ^ + 2 0 O (0 2 + 20,) - 40, 2 . 


The sign of the square root is determined from the initial conditions of the problem. X itself is cer- 
tainly positive if the motion of the satellite is direct, for, by definition, 


x = 0 O + 0 2 - 6k<f 2 - r 2 + 2(^T)-7jf) . 


On the other hand, if, in the fixed xy-system, xy-yx = r 2 0is the velocity at a given instant, 


~ r 2 (0- 1) , 
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so that 


k = r 2 ( 2^ + 1 ) . 

Therefore k is positive even for > -1/2 . 

The advantage of the differential equation (Equation 20) is that it contains only one variable, so 
that the expression for the perturbation (on the right side) also depends only upon r, f, r and c. 

In the application of the methods of numerical integration, it suffices, therefore, to set up a single dif- 
ference table. This step could be of particular importance when working with electronic computers. 
The relatively complicated structure of the disturbing function does not constitute an appreciable 
barrier for the computers. 

It must also be mentioned that Equation 20 gives the function r(t) independently of the coordi- 
nate system selected, and that therefore its solution is as valid in a fixed as in a revolving system. It 
is here that this solution differs from that which Hill gave for the system*. Hill's solution was lim- 
ited to simple-periodic orbits in a revolving t) system. Although the Hill experiments yield — of 
the known periodic inequalities — only the so-called "variation" which is periodic in the revolving 
system, a similar analysis of the differential equation (Equation 20) would also yield the "large in- 
equality" of the moon motion (largest periodic term of elliptical motion), the "evection," and some 
other features. Likewise, considering the secular perturbations, we find the apsidal motion, but it 
is inherent in the problem that the terms due to the solar parallax, the eccentricity and inclination 
of the earth's orbit do not appear. 

Provided that we have succeeded in determining r = r ( t ) — either by numerical integration of 
Equation 20 or by a theoretical formula (developed perhaps through a sufficiently converging series)- 
it is always possible to determine the corresponding coordinate 4>[t) by means of a simple quadrature. 
By letting <f = r cos 4> and 77 = r sin <£, Hill’s equations (Equation 1) take the form: 

(r “ r<p 2 - 2r <p) cos <p ” ( 2r<p + r<p + 2r ) s in <p ~ ^3 - r cos <p , 

(r “ r<p 2 ” 2r<^)sin0 + (2f <p + r <p + 2f ) cos <p ~ rsin<£ 


from which 


r “ rep 2 - 2r<p = 3r cos 2 <p r , 

r * 

. .. 3 

2r<p + r<p + 2r = “ 3r s in cos <p - — ^ r s in 2<£ 




> 


( 21 ) 


•Hill, G. W., "Researches in the Lunar Theory,* in The Collected Works of George William Hill, V ol. 1, Washington, D. C.: Carnegie 
Institution, 1905* 


7 



follow. Hill has shown from the second equation of Equation 21 that r = r (t ) is found by quadrature 
if 0(t)is known, for it follows from 

2 —(<j> + 1 ) = -Lp + ~2 s in 2<fi\ 


that, with k constant, 


n 


0 + i 


3 

exp - 


\ sin 20 

JTTT dt 


In the same manner we also find from the first equation of Equation 21 that 


0 2 + 20 


~F + ~ - 3 cos 2 0 


or 


0 + 1 


/TTT 

r r r 3 


1-3 cos 2 0 


( 22 ) 


from which it follows, by integration, that 


t - t. 


r 


d0 


i/TTX 

r _ 3 


+ 1-3 cos 2 0-1 


(23) 


The above relation yields t = H<p) and, therefore 4> - <£(t). 

If 4 > = < t> 0 + nt = 4> 0 + 1 is the true length of the satellite in a fixed coordinate system, then 

* = 4 > + 1 = 4 . 


for undisturbed motion, if p denotes the parameter of the orbit ellipse; we also have 


therefore 



4> + l 



( 24 ) 
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This differential equation can therefore be considered as an approximation for Equation 22 as long 
as the term (1-3 cos 2 <p) in the radicand of the square root remains small by comparison with the 
other terms. For the earth's moon, which is of particular interest here, this condition is still satis- 
fied, for, since the unit of length is so selected that the mean earth- sun distance is about 69, the 
average value of r (normally 390 times smaller) is about 3/175 the expression { r/r ) + (l/r 3 ) is there- 
fore of the order p/r 4 or l/r 3 , or about 180, whereas the attachment term varies from -2 to +1. 


INTRODUCTION OF A NEW INDEPENDENT VARIABLE 

The form of the differential equation (Equation 20) suggests a solving process which is a 
modification of that used in the undisturbed problem (Equation 18). If a new variable q is in- 
troduced instead of time t by means of q = l/r and the initial condition q(t 0 ) - 0 so that q disap- 
pears simultaneously and increases monotonically with (t - t 0 ) , we find that* 



(25) 





With the above equations, the undisturbed equation of motion (Equation 18) assumes the form 


1 


r 


2 


t 



= 0 . 


But since 


d 

dq 




(26) 


(27) 


is constant, and Equation 26 goes over into the integrable form: 

r 1 " + a 2 r' = 0 . 


(28) 


•The derivatives of q are denoted by primes. 
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On the other hand, in the disturbed problem (Equation 20) we have 


r 


m 


+ 


1 


rt 


r 


r 


4>(q) = <£(r , r f , r",C) 


or 


r w + a 2 r ' = <#>(q) , (29) 

where ^(q) represents the expression for the perturbation (right side of Equation 20) as a function of 
q. And because of Equation 27, 



4>(g) 


or 


f 1 

Jo r 


, 

(q) d q • 


(30) 


Therefore, if we set 




<fi + r' 



(31) 


we get 


t "' + a 0 2 r' - g(q) 


(32) 


where - (l-r 0 ")/r is a constant which is known from the initial conditions. 


If g(q) is a sufficiently small perturbation in the interval |qj <o, the integration method de- 
scribed in a companion paper* could be utilized. If, for s = a Q q, we set c 0 (s 2 ) - cos s, 
Cj(s 2 ) = (sin s )/s, c 2 (s 2 ) = (l - cos s )/s 2 and c 3 ( s 2 ) = (s - sin s)/s 3 , •••, then the recurrence 

formula 


( 33 ) 


"Stumpff, K., "Calculation of Ephemerides from Initial Values," NASA Technical Note 0-1415, 1962. 
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and the differential equation 


dq 


( c „ + i q n+1 ) = c n q n 


(34) 


are applicable. Therefore the Taylor development 


rfq) = r 0 + yr r 0 ’ q +-^T r"q 2 + jrr^'q 3 + 


can, by replacing the reciprocal factorials with the expressions given by Equation 33, be written 
rfq) = r 0 + u ( c i + s * c j)<> + r o" ( c 2 + s2c 4 • 


If we again set s 2 = a 2 q 2 and arrange according to powers of q, we have 


r(q) = r 0 + r 


o' C,q + r 0 "c 2 q 2 + (r“ + aj r 0 ') c 3 q 3 + a„ 2 r 0 ")c 4 q 4 


a series which, because of Equation 32, becomes 


rfq) = r„ + r 0 'c,q + r 0 "c 2 q 2 + g„c 3 q 3 + g^q 4 + g 0 " c s q 


5 + 


(35) 


Now we set 


r 0 q = zt , (36) 

where r = k(t -t 0 ) denotes the intermediate time expressed in units of 1A days; we introduce the 
expressions y lf y 3 ? (which are invariable against coordinate transformations) through 


l-n = r n 2 r, 


y*' 


and we set 


y \ T V\ * y 2 T *" ^ 2 * y v? ' 


Equation 35 now assumes the form 

r(q) = r 0 [l + c^jz + c 2 75 2 z 2 + c 3 77 3 z 3 + •••] , (37) 

a series, the terms of which are small above the third order, and which, in practical cases and for 
moderate intermediate times, converges very rapidly. The association between r and q is provided 
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by the differential equation 


dr 


r dq 



+ c , V 1 Z + C 2 T)j Z 1 + 



But 


r oJ C n % Z " d( l ~ r 0^nJ C n (zT) " df l = r o" + ‘ f C n <T 


dq 


and therefore, according to Equation 34, 


(38) 


» 


-j 


c 77 z n dq 
n 'n 1 


r 


n + l 
0 


c 

n 


+ i q 


n + l 


y n c n 


T v n 


If we divide the integral of Equation 38 by r, we obtain the main equation 

1 = z + c lVl z 2 + c 3 t, jZ 3 + c 4 7j 3 z" + . ( 39 ) 

For undisturbed motion (g = 0), it limits itself to the first three terms, since v 3 , v 4 , ■■■ disappear. 
The c -functions 


C n(^ 2 ) 


J_ Xi_ x 4 

n! " { n + 2 ) ! + (n+4)! _ 


have the always- real argument 


(vO J 


or, since 



CLOSING COMMENTS AND OUTLOOK 

The aforementioned integration method replaces the numerical integration method of Hill’s dif- 
ferential equation (Equation 1) for the case where we can, as a first approximation, consider the 
satellite motion as a Keplerian ellipse with any desired eccentricity. This method permits the 
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determination of r(t) directly through iteration - at least under some conditions, for a certain fairly 
extended time interval in the vicinity of the initial time t 0 . The iteration process itself is limited to 
the solution of the transcendental main equation (Equation 39), since the quantities g 0 , g 0 ', • • • can be 
derived from the initial conditions. 

Here we have attempted to demonstrate this method on a problem — Hill’s moon problem — which 
perhaps has little practical importance, but is simple enough that its solution may be written explic- 
itly and understandably. Thus it serves as an example for other, more difficult problems, which can 
be handled with the same general principles and solved in the same manner. This method can be 
used to attack certain special cases of the probleme restreint, in which conic-section orbits can be 
considered good approximate solutions. For example, we could consider the case of the sufficiently 
close orbit of an infinitesimal body around one of two finite masses; i.e., either the actual satellite 
problem in which the massless body moves around the lesser mass (planet) or the problem of dis- 
turbed planetoid orbits, when the motion occurs about the greater mass (sun) and the distance from 
the disturbing planet remains sufficiently large. 

In these cases it is also possible to apply the aforementioned method rationally, since we are 
again concerned with a problem of motion of the fourth order, which can be reduced to the third 
order by applying the Jacobian integral. But here the difficulties are considerably greater, although 
not insurmountable. The elimination process described becomes much more complicated because 
the quantities to be eliminated (<f, T: , k , y) appear not only in the distance and velocity of the 
body in reference to its central mass, but also in the distance and velocity in reference to the dis- 
turbing mass. In Hill’s problem, the last step of this process leads to the elimination of X from the 
two algebraic equations (Equations 13 and Equation 14) which are of the second and fourth degrees. 
The resultant equation (Equation 15), which solves the problem, is therefore obtained in the form of 
a six- row Sylvester determinant which is set equal to zero. The same process applied to the prob- 
leme restreint leads to a correspondingly-to-be-determined quantity k that is to be eliminated from 
the two algebraic equations of the 10th and 14th degree, so that the final solution appears as a 24- row 
Sylvester determinant which is equated to zero, the strict solution of which would be hopelessly com- 
plicated. However, as long as we restrict ourselves to those cases in which the motion of the satel- 
lite or planetoid can be viewed as an only slightly disturbed Kepler orbit, the elimination can be 
accomplished with the help of a rapidly converging approximation process. The demonstration of 
this must be reserved for another treatise. 


NASA -Langley, 1963 
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